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ABSTRACT 

We determine the dependence of the initial helium abundance and the present- 
day helium abundance in the convective envelope of solar models (Y ini and Y sur f 
respectively) on the parameters that are used to construct the models. We do so 
by using reference standard solar models to compute the power-law coefficients of 
the dependence of Y ini and 3^ ur f on the input parameters. We use these dependen- 
cies to determine the correlation between Y ini and F sur f and use this correlation to 
eliminate uncertainties in Y- m i from all solar model input parameters except the 
microscopic diffusion rate. We find an expression for Y- m i that depends only on 
y sur f and the diffusion rate. By adopting the helioseismic determination of solar 
surface helium abundance, Y| urf = 0.2485 ± 0.0035, and an uncertainty of 20% 
for the diffusion rate, we find that the initial solar helium abundance, Y™ 1 , is 
0.278 ±0.006 independently of the reference standard solar models (and particu- 
larly on adopted the solar abundances) used in the derivation of the correlation 
between Yj n i and Y sm f. When non-standard solar models with extra mixing are 
used, then we derive Y"™ 1 = 0.273 ± 0.006. In both cases, the derived Y™ 1 value 
is higher than that directly derived from solar model calibrations when the low 
metalicity solar abundances (e.g. by Asplund et al.) are adopted in the models. 
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Subject headings: Sun: abundances — diffusion — Sun: helioseismology — Sun: 
interior 



INTRODUCTION 



The primordial helium abundance, Yp, can be determined to a high accuracy thanks 
to the measure ment of cosmological parameters, particularly of the bar yon-to-photon ratio 
7/10, by WMAP JSpergel et alil2007h . The most recent WMAP res ults Larson et al.l lioioh 
give Yp = 0.2486 ± 0.0006 (using the analytic formula presented by ISteigmanl 120071 relating 
Yp and rjio), but rely on the assumption of Standard Big Bang Nucleosynthesis (SBBN). 
On the other hand, empirical determinations of Yp, that can be used to test the SBBN 
scenario, rely on the determination of helium abundances in metal poor HII regions and 
extrapolations to zero metalicity, assuming a linear relation betwe en helium and metalicity 
(or oxygen) with a given slope AY/AZ (see IPeimbert et al.l 120071 for a detailed account of 
the meth od and uncertainties i nvolved). The results are not exempt of controversy. For 
example, IPeimbert et al.l (120071 ) find Yp = 0.24 77 ± 0.0029, in excell ent agreement with 



WMAP+SBBN results, while other authors (e.g.. Ilzotov &: Thuanll2010l ). find higher values, 
Yp = 0.2565±0.0010 (stat)±0.0050 (syst), and suggest deviations from SBBN. Discrepancies 
are not confined to Yp, but are also found in the given values of AY/AZ. Other attempts to 



determine AY/AZ, both observationall y and theoretically, inc 



ing of the lower main sequence (see e.g. ICasagrande et al.ll2007 andlGennaro et al.ll2010l for 



naries (Ribas et al. 



recent works), the study of galactic stellar populations ( Renzinil 



ude the use of the broaden- 



19941 ). modeli ng eclipsing bi 



200 0)), construction of galactic ch emical evolution models (ITimmes et al. 



19951 ; IChiappini et al.ll2003l : ICarigi fc Peimbertll2008l ) . The values obtained for AY/ AZ range 
approximately between 1.4 and 2.5 for different authors and methods. 

The initial composition of stars determines their structure and evolution and, by ex- 
tension, the characteristics of stellar populations. The initial helium abundance of stars 
is in most cases impossible to determine directly, hence indirect methods have to be used. 
A common approach used in stellar and stellar populations modeling is to assume a value 
for AY/AZ and use the solar helium and metalicity as an anchor point for the relation. 
This is a routine procedure used in, for example, libraries of stellar evolutionary tracks to 



defin e the initial composition used in the stellar models (e.g., ISchaller et al.lll992t lYi et al. 



200ll : iPietrinferni et al.l 120041 ) The initial solar helium abundance, however, can not be de- 
termined directly. It is well established that element diffusion takes place in the Sun and 
that the present-day photospheric abundances of helium and metals are lower than the 
initial ones. Therefore, the initial solar helium abundance is usually obtained from cal- 
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ibrating solar models, which then introduces model dependencies on the determined val- 
ues. An illustrative example is the effect of solar photospheric metalicity determinations 
on the inferred initi al solar helium abundan ce. Solar models that use o lder solar abun- 
dances like those of iGrevesse Noelsl (119931 ) or iGrevesse & Sauvall (119981 ) l ead to models 
w ith initial helium values in the m odels in the range Y m i ~ 0.272 — 0.275 (IBahcall et al. 



200 ll ; iBoothrovd fc Sackmann 



200 31), while models that incorporate newer deter minations 



of solar metalicity (lAsplund et al. 
Y ini « 0.260 - 0.263 feahcall et al. 



2005 



2005 



adders et a. 



Guzik et al. 



2009 



2005; 



Asplund et a 



Serenelli et al. 



2009 



2009 



e.g.) give 
While the 

models with the older abundances have present-day surface helium abundances F sur f consis- 
tent with helioseismic determinations of the helioseismology determinations of the present- 
day solar convection-zone helium abundance, the models with the lower metal abundances 
do not. It is interesting to note also that if we adopt Y ini w 0.260 — 0.263 values together 
with initial solar metalicities derived from the same solar models, i.e. Z ini « 0.014 — 0.015, 
then values as low as 0.8 are obtained for AY/AZ (assuming Yp from SBBN). This is in 
disagreement with all accepted determinations of AY/AZ as we have summarized above. 

In this work, we propose a simple method to determine the initial solar helium abun- 
dance based on the use of solar models to calibrate the relation between the initial and the 
present-day surface helium abundances. Although derived with the help of solar models, we 
then show that the predicted initial solar helium abundance (Y™; not to be confused by the 
initial helium abundance Y ini resulting from solar model calibrations) is a robust result, i.e. 
a very wide range of solar models (standard and non-standard) lead to the same Y™ value. 
The structure of the paper is as follows: in §|2]we determine the effect of changes in the input 
parameters used in solar models (including metal abundances, diffusion rate, nuclear cross 
sections, etc.) on the initial and present-day helium abundances of the solar models; in § [3] 
we derive a simple analytic relation between Y [ni and y sur f that minimizes the uncertainties 
originating in the inputs to solar model calculations; in § H] we present our main results and 
discuss their range of applicability; finally, we summarize our findings and conclusions in § O 



2. POWER-LAW DEPENDENCE OF F surf AND F ir 



For this pa per, we use t he same technique that has been used earlier for analyzing 
neutrino fluxes (jBahcalllll989l see Chap. 7). We express the dependence of the initial (F ini ) 
and present-day (F surf ) helium abundances of solar models as a product of powers of the 
different input parameters: 
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In these equations, the input parameters in solar models (nuclear cross sections, abun- 
dances of individual elements, age, etc.) are represented by a set of parameters {pi} i=1 N with 
central, or reference, values denoted with the sub-index "0". These parameters comprise: 
metal abundances, nuclear reaction rates, element diffusion rate, solar surface parameters 
(luminosity, age); a detailed list is given in Table [TJ The exponents and (3i respectively 
determine the dependence of Y ini and F sur f on the input parameters. They are defined as 

<91ogF ini d\ogY sur{ 

Ui = -7T] , Pi = -TT, • (3) 

o log pi a log pi 



The indices in Equations (J3J) were calculated for standard solar models (SSM) using 
different determina t ions o f pho tospheric solar compo sition. We chose the compilations by 
Grevesse fc Sauvall (119981 ) and lAsplund et al.l (120051 ). hereafter GS98 and AGS05 r espec 



tively, and conservative uncertainties, according to the definition by iBahcall et al.l ( 12006 
hereafter BSB06). A discussion of the central values for all parameters, as well as their 
uncertainties, is presented in that work. To derive the indices we proceeded as follows: for 
each input parameter pi we computed a set of 50 SSMs where the value of pi is drawn from 
a Gaussian distribution, with standard deviation equal to the assumed uncertainty of the 
parameter, while all other input parameters are kept fixed and equal to their central value. 
The logarithmic derivatives defined in Equations [3] are then obtained from fits to each of the 
small Monte Carlo simulations described above. Some examples are shown in Figure HJ The 
values of these dependencies are tabulated in Table [H The exponents derived for both GS98 
and AGS05 solar models are similar. This is so because relative changes in input parameters 
produce relative changes in solar model properties that do not, to the first order, depend on 
the reference solar model. Even more important for this work is that the ratios are 
almost independent of the reference solar models used. Differences in the values obtained for 
GS98 and AGS05 models will be used later as a way to determine the systematic uncertainty 
of the method presented in this work. 

Radiative opacities can not be represented by a single parameter. Hence to account 
for the influence of opacities on Y m i and y sur f we have proceeded in the following way: for 
both AGS05 and GS98 compositions we have computed two solar models, one with OPAL 
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opacities and the other one with OP opacities, while keeping all other input physics the 
same in both models. For each composition we define the la uncertainties in Yj ni and K, ur f as 
half the difference in the values obtained for models with the same composition but the two 
different opacity sources. We find that for F sur f the uncertainty from opacity is 0.23%(0.20%), 
while for Yj ni is 0.30%(0.28%) for the AGS05(GS98) composition. The uncertainties between 
Yi n i and Ygurf are found to be fully correlated. We account for opacity uncertainties by adding 
Gaussian dispersions with the above defined standard deviations to the distributions of Y ini 
and Y suri that we generate using Equations [T] and [2J 

We have verified that Equations [1] and [2] describe the behavior of Y{ D i and Y sur f in solar 
models very well by comparing distributions generated with these equations against results 
obtained with full solar models. In particular, we have used the large set of solar models 
computed for the Monte Carlo simulation done by BSB06, where different assumptions for 
the reference solar composition and abundance uncertainties were adopted. Results of the 
comparison are presented in Figure [2j The panel on the left shows results for the GS98 
composition with conservative (large) abundance uncertainties (see BSB06 for details); while 
the panel on the right shows results for models using AGS05 composition and much smaller 
uncertainties, as originally given by AGS05. In both cases, results from the full solar models 
obtained from the BSB06 Monte Carlo simulations and the resulting distributions for the 
Yini and Y" sur f obtained by means of Equations [1] and [2] are shown. We see in both cases that 
results from full solar models are reproduced very closely by the power-law expressions. 

Having verified the validity of power-law expressions in reproducing the correct distribu- 
tions of Yi n i and Y" sur f , in what follows, we concentrate on using these expressions to constrain 
the solar initial abundance, Y™ 1 , in a model independent manner. We do this by using the 
helioseismically determined solar convection-zone helium abundance. 



3. CORRELATION BETWEEN F surf AND Y ini 

In this paper we assume, unless specified otherwise, that only microscopic diffusion of 
helium and heavy elements modify the composition of the outer layers of the Sue0. Therefore, 
in the absence of diffusion, Y" sur f would exactly match Yj ni . The most evident effect of diffusion, 
through gravitational settling, is to decrease Y" sur f with respect to Y ini . Thus Ys Ur f and the 
diffusion rate are anticorrelated. Diffusion also modifies the composition in the solar core. 
In particular, gravitational settling of metals in the core leads to a higher radiative opacity 



lr The surface abundance of lithium and beryllium are also modified by nuclear burning but that is not 
relevant given the scope of the present work. 
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and hence, to a steeper temperature profile. The slightly increased temperature in turn 
increases the rates of nuclear energy production. However, the solar luminosity is fixed. Thus, 
the increase in energy production due to the larger temperatures has to be compensated 
by decreasing the amount of available hydrogen or, equivalently, by increasing the helium 
abundance in the core. As a consequence, Y m i is positively correlated with the diffusion 
rate. We call these two combined effects, the negative correlation of Y SUT f the the positive 
correlation of Y ini with the diffusion rate, the 'direct influence' of diffusion. 

As can be seen from the coefficients listed in Table [TJ the direct influence of diffusion is 
the only effect that modify 3^ ur f and Yj n i in opposite ways. For this reason, in what follows 
we treat the effect of diffusion separately. Thus Equations. [1] and [2] are modified to: 
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n 

i=l,AT;i^Diff 

n 

.i=l,JV;i^Difi 



Pi 
Pifi 

Pi 
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(4) 



(5) 



Using the expressions in brackets in the equations above, we obtained distributions for 
K; Ur f and Y ini that exclude the direct effect of diffusion. The correlation between these two 
quantities is now very tight as can be seen in Figure |3j Naively, we would expect no dispersion 
to be present, since the diffusion rates are kept fixed. However, keeping the diffusion rates 
fixed does not imply the total amount of helium (and metals) diffused out of the convective 
envelope is constant. The effective rate of diffusion in the solar envelope is directly related 
to the extension of the convective envelope which acts as a reservoir of helium. Therefore, 
any changes to the model inputs that have an effect on the mass contained in the convective 
envelope will be reflected in a deviation of the behavior of Y" sur f with respect to that of Y ini 
(i.e., effectively changing the values of the power-law exponents «j and 0i to a small extent), 
leading to the dispersion seen in Figure |3j 



Following a procedure analogous to that used by lHaxton fc Serenellil (120081 ) to study 
solar neutrino fluxes, we take advantage of the small dispersion in the correlation between 
y sur f and Yi n i to express the uncertainties in Y ini in terms of those in Y Bm { while minimizing 
the uncertainty in all parameters except for diffusion. To do this we combine Equations H] 
and |5] to obtain: 



± ii 
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Pi 
Pifi 
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where 



7i = a* - if(i,s)A 



(7) 



The constant ifri s) is obtained by doing a linear fit to the simulated data for log(Yjni) and 
log(Kj Ur f) using the power-law expressions while ignoring the direct dependence on diffusion. 
For the GS98-C case we obtain i^(i,s) = 0.836. The residuals of the fit have a dispersion of 
0.4% (left panels of Figure[3]). For the exponent in the diffusion factor we obtain 7 Diff = 0.088. 
We can now write 



Y ■ 



surf 



Y 



surf ,0 



0.836 



gDiff 
PDiff.O 



[1 ± 0.4%(ND)] 



where the uncertainty of 0.4% represents the total uncertainty due to all non-diffusion param- 
eters (note that only abundance- uncertainties play an appreciable role). If we used optimistic 
abundance-uncertainties instead, the total non-diffusion uncertainty would be only 0.2%. 

To proceed further, we need to consider the uncertainty of the diffusion rate (hereafter 
CDiff)- The best agreement between the structure of solar models and the helioseismically 
determined sol ar structure is obtained when helium and heavy-elemen t diffusion is included 
in the models (iBahcall et al.lll995l : IChristensen-Dalsgaard et al.lll996l ). However, the actual 
diffusion rates of helium and metals a nd, in particular, the associated uncertainties are not 
well constrained. iThoul et al.l (I1994J ) have esti mated an uncert a inty o f about 15% in the 
diffusion rates. A similar result was obtained by iFiorentini et al.l (119991 ) by using conditions 
at the base of the co nvective zone to constrain the diffusion rates. A detailed study by 
Turcotte et al.l (119981 ) shows that diffusion velocities depend on the different assumptions 
made in the calculation for different elements and, at different regions of the Sun, the vari- 
ation ranges from a few percent to about 40%. Here, we adopt an intermediate fiducial 
value (TDiff = 20% as the uncertainty in the diffusion rates. It should be noted that the 
uncertainty in Y^ because of uncertainties in the diffusion rate are only one-twelfth that of 
uncertainties in the diffusion rate. Then, for o"Diff = 20%, the contribution of diffusion to 
the total uncertainty in Yj ni is only 1.7%. Thus, 



Y ■ 

1 mi 

Yni,0 



Y s 



surf 



Y 



surf ,0 



0.836 



[l±1.7%(Diff)±0.4%(ND)] 



(9) 



We now scale the relation to the Sun by adopting the present-day surfac e helium abun 
dance determined using helioseismology. We adopt the value determined by iBasu fc Antia 



(12004J ). i.e., Yq = 0.2485 ± 0.0034, where the uncertainty a ccounts for svstemat ics, includ- 



ing those caused by uncertainties in the equation of state (see lBasu fc Antiall2008l for a recent 



discussion about the determination of Y"q ). Then, the solar initial helium abundance Y c 



mi 
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can be expressed as 

/O 2485 \ a836 

Y™ = >ini,o [1 ± l.l%(Helio) ± 1.7%(Diff) ± 0.4%(ND)] . (10) 

In the above expression, standard solar models (SSMs) play a two-folded role. First, a 
reference model is used to get the scaling factors Yi n i,o and Y^ U rf,o- Second, models are used 
to determine the exponent that relate Y"Q Urf and Kjurf,o- However, as we show in the following 
section, Equation [10] has predictive power independent of the (standard) solar model used 
as a reference. 



4. RESULTS 

As a fir st test of Equa t ion ITUl we use values 3^ n i,o and 3^ U rf,o obtained from SSM cal- 
culations of ISerenelli et al.l (120091 ) . The results are summarized in Table [2] where we have 
identified the solar models by the solar composition adopted for each of them. The second 
and third columns give the SSM predictions for the initial and present-day surface helium 
mass fraction. The fourth column gives the Y"g m values estimated using Equation [TO) the 
uncertainty in all CcLSGS IS (Xy-ini = ±0.006. Comparing results for the different standard solar 
models we find that the dispersion in Y^ 111 is about ten times smaller than cry-mi and hence 
we conclude that the estimated values of are consistent with each other, and they are 
independent of the standard solar model used as reference. 

If standard solar models based on AGS05 composition are used to derive the relation 
between Y" ini and Y sm f, we get Kng\ = 0.903 instead of 0.836. We have used GS98 and 
AGS05 for this comparison in an attempt to bracket what the real solar abundance probably 
is. These two compilations represent the two ends of the currently accepted range of values 
for solar photospheric abundances q Therefore, we also expect that our estimates of Kn t s) 
for the GS98-C and AGS05-Op cases represent the two extremes of the possible range of 
values for this qu antity. Replacing Kn ^ = 0.903 in Eq. [10] and again using the standard 



solar models from ISerenelli et al.l ( 120091 ). we get the YZ 111 values listed in the last column of 



Table [5] A comparison of the Y^ 1 results obtained with the two values of i^(i,s) for each 
SSM shows that the uncertainty in the determination of if(i,s) has a very minor effect in the 



2 In fact, newer determinations of solar abundances, e.g. by lAsplund et al.l (|2009l ); ICaffau et all (|2010l ) . 
point towards abundances o f indiv idual elements in the solar photosphere that are systematically higher 

J2005|) (although consistent within the quoted uncertainties), but lower than 



than those in 



Asplund et al 



those in iGrevesse fc Sauval 



199 
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predicted Y"g m values. The largest difference is found for the AGSS09ph model, for which 
AYg m = 0.0017, a factor of ~ 4 smaller than <7y™ estimated above. 

We have shown that the predicted Y™ values have a very small dependence on the solar 
model used to determine it. The choice of standard solar models for deriving the exponent 
Kng) also introduces a small scatter in the resulting Y™ 1 values. Therefore, we average all 
values Yq 1 presented in Table [2] to determine the final result and also use the scatter to 
define the systematic uncertainty of the method. The final result that we obtain is 

Y™ = 0.278 ± 0.006(helio, param) ± 0.002(syst) (11) 

where the first part of the uncertainty includes contributions from helioseismology and solar 
input parameters, and the second one is the systematic uncertainty of the method as defined 
above. 

We have tested the robustness of this result by applying Equation (fTOj) to a variety of 
solar models, standard and non-standard, i.e. including additional physics not considered 
in SSMs. As a first consistency check, we have applied Equation [10] to the two large sets of 
SSMs computed by BSB06 used to assess the validity of the power-law expansions obtained 
in §[2j Averaging all results for the AGS05-Op Monte Carlo set presented in BSB06, we find 
Y™ = 0.2780 ± 0.0043 and for the GS98-Cons set we get Y™ = 0.2783 ± 0.0041, i.e. results 
fully consistent with the basic derivation of Y^ 1 and its uncertainty based on the four SSMs 
presented in Table [2j These uncertainties do not include errors in the helioseismic results. 

To further test the robustness of our derivation of Y™ we have resorted to results for 
standard and non-standard solar models from the literature. Non-standard models have 
been considered only when the modification to the input physics results in a model that is 
in better agreement with helioseismic results than its SSM counterpart. As a last test, we 
have computed additional non-standard models to test the influence of turbulent diffusion 
on the robustness of our predictions. 

Here is a summary of our observations: 



Turcotte et al.l (119981 ) have tested the effect of different diffusion schemes (including 
radiative levitation in some cases) and the use of monochromatic opacities. By applying 
Equation (fTOl) to their models that include diffusion, we obtain an average Y™ = 
0.2780 ± 0.0017, where the uncertainty refers to the scatter between models. From all 
their models we infer Y™ 1 consistent within la of our result in Equation [TTJ 



- Some standard and non-standard models have been presented in lBahcall et al.l ( 1200 ll ). 
From these models (except their models with no diffusion and with a mixed core; 
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both clearly ruled out by helioseismology) we determine an average value 3^ ni 
0.2774 ± 0.0018 (scatter between models). The model that most deviates from this 
results includes rotation^ which, by introducing additional mixing, tends to inhibit the 
efficiency of diffusion. For this particular model we get Y™ 1 = 0.2725, still within la 
of the predicted value in Equation [TT1 



Boothroyd fc Sackmannl (120031 ) have computed a large set of different solar models with 
different variations in the input physics and parameters. For all these models we derive 
Yq 1 consistent within la of the central result shown in Eq. [TTJ The only exception is 
their model where the diffusion velocity of helium is increased by 20%, for which we 
derive FQ Urf = 0.2873, a 1.5a difference with our derived central value for Y^ urf ; but 
this model shows a convective envelope too shallow compared to the helioseismic value, 
and hence, we do not consider the discrepancy to be important. Averaging all results, 
we get Y™ = 0.2783 ± 0.0013. 

We have computed solar models th at include turbulent diffusion as given by Equa- 
tion (2) of lProffitt fc Michaudl ( 199ll ). where the diffusion coefficient is D T = D [pcz/p] 3 , 
and p is the dens ity and pc,7, the density at the base of t he convective envelope. 
As pointed out by IChristensen-Dalsgaard fc Di Maurol (120071 ) , even a small value for 
D (200 cm 2 s _1 ) eliminates the bump in the sound-speed difference between mod- 
els and the Sun below the convective envelope. We computed models for D = 
300, 1000, and, 3000 cmV 1 with both GS98 and AGS05 compositions. Applying 
Equation flTD]) to these models, we find = 0.2725, 0.2690, 0.2665 for increasing D 
and irrespective of the composition adopted. It should be noted that all models with 
turbulent diffusion predict sound-speed profiles with average root-mean-square devia- 
tions from helioseimic results that are comparable to the rms deviations of SSMs. In 
fact, for models based on the GS98 composition we get an average rms sound speed de- 
viation of 0.0009 for the SSM while for models with D = 300, 1000, and, 3000 cmV 1 
we get 0.0009, 0.0009, and 0.0010 respectively. For models based on the AGS05 com- 
position, results are 0.0047, 0.0053, 0.0054, and 0.0055 for the SSM and models with 
increasing D as given above. This is a result of the fact that, although for models 
with turbulent diffusion the bump below the CZ disappears, the lower central helium 
abundance in the models lead to a slight worsening of the agreement with the solar 
sound speed towards the center. This is also reflected in the frequency separation ra- 
tios, particularly r 02 (IBasu et al.l 120071 ). which in all three cases show worse agreement 
with the observed solar r 02 than the corresponding SSM constructed with the same 



3 According to iPinsonneault et al.l (|1999i ). this model sets a reasonable upper limit to rotation- induced 
mixing, compatible with the observed lithium depletion in the Sun. 
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composition. Larger effects in this direction are found for increasing Dq values. From 
these results, we conclude that non-SSMs accounting for the effect of extra mixing lead 
to the prediction: 

= 0.273 ± 0.006(helio,param) ± 0.002(syst), (12) 
consistent with the result from SSMs given in [12j 

5. SUMMARY AND CONCLUSIONS 

We have used standard solar models to study the dependence of the initial (Yni) and 
surface (Y sur f ) helium abundances on different input physics and element abundances entering 
solar model calculations. As expected, there is a very tight correlation between these two 
quantities, particularly if direct effects of diffusion are not taken into account. We have 
used this tight correlation to express the initial solar helium abundance (Yq 111 ; not to be 
confused with the initial value predicted by models) as a function of the solar surface helium 
abundance, YQ mf , determined from helioseismology, while minimizing the uncertainties from 
solar model input parameters. This result is expressed in a compact form in Equation [TOJ 
A slightly more general expression where the scaling of the helium and metals diffusion rate 
is still explicit is given in Equation |HJ 

We have applied Equation (fTUj) to different, recently published standard solar models 
listed in Tabled and derive Y"™ 1 = 0.278±0.006 (Equation [11]). The uncertainty is dominated 
by uncertainties in the diffusion rate of elements and those in the helioseismic determination 
of the solar surface helium abundance. The dispersion in the results for the different solar 
models is one order of magnitude smaller. The results are also robust in terms of the standard 
solar models used to derive Equation (TDJ 

To test further the robustness of our results, we have compiled results from literature 
for standard and non-standard solar models (though only those non-standard solar models 
that improve agreement with helioseismic inferences of solar structure). In all cases, even for 
the non-standard ones, applying Equation [TU] to results from solar models yield Y™ results 
consistent to within la of our central value of Y™ 1 = 0.278 ± 0.006. 

Finally, we have also computed solar models including turbulent diffusion as a phe- 
nomenological approach to eliminate the bump in the sound-speed difference between mod- 
els and the Sun that is seen below the convective envelope. For these models, a small value 
of the free parameter Dq, Dq = 300cm 2 s _1 , leads to Y"g m = 0.2725, consistent within la 
with our prediction given in Equation [TTJ Larger values of Do predict systematically lower 
values for Y™ 1 but lead to slightly degraded agreement with helioseismology, as evinced by 
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frequency separation ratios (probably a result of the lower central helium abundance pre- 
dicted by these models). Thus, we find no evidence, at least from helioseismology data, that 
support the need of larger D values. It has to be kept in mind that this formulation of tur- 
bulent diffusion represents only a phenomenological approach to eliminating the bump in the 
sound-speed difference profile. More realistic mode l s that account for rotationally-induced 
mixing, like that one derived by iPinsonneault et al.l ( 119991 ). lead again to Y™ 1 = 0.2725. 



We conclude that the initial solar helium abundance, as inferred from the present-day 
solar surface helium abundance determined by helioseismology, is = 0.278 ± 0.006. Al- 
though solar models are used to reach this conclusion, we find this result is almost indepen- 
dent of which solar models are used in its derivation. If Y™ 1 is determined from non-standard 



solar models that account for the effects of extra mixing, we obtain Y™ 1 = 0.273 ± 0.006. In 
all cases, i.e., for both standard and non-standard models, we infer Y™ 1 values that are higher 
than the initial helium abundance obtained in solar models that use the solar abundances 



by lAsplund et al.l (120051 ) or the most recent determination by lAsplund et al.l (120091 ) . This 
result points towards a deficit in solar models using these abundances; whether the problem 
lies in the abundance determinations or in the constitutive physics of the models is beyond 
the scope of this work. 

Depending on the solar metalicity adopted, we derive a helium-to-metal enrichment of 
AY/ AZ ~ 1.7 — 2.2, in line with standard predictions of Galactic chemical evolution and 
other determinations of AY j AZ available in the literature. 



This work is partly supported by NSF grant ATM-0348837 to SB. 
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Fig. 1. — Some examples of the derivation of power- law exponents. Left panels show results 
for y ini and right panels for Y SVLr {. From top to bottom, results correspond to changes in the 
diffusion rate, carbon, and iron abundances. Black filled circles: results for GS98 compo- 
sition; red empty squares: results for AGS05 composition. In each panel, lines denote the 
linear fits to the model results; the slope is the power-law exponent (see Equations EJ). 
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Fig. 2. — Comparison between Y sur { vs. Y m i relation obtained with full solar models (big 
black dots) with the obtained using power- law expressions given in Eqs. [1] and [2] and Tabled] 
(small red dots). Results for GS98 composition and conservative uncertainties and AGS05 
composition with optimistic uncertainties are shown. 
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Fig. 3. — Y ini vs. F surf distributions obtained using power-law expansions (Equations 0] and 
but neglecting variations in the diffusion rates (first factor in the equations). 



- 18 - 



Table 1. Power-law exponents for solar model input parameters. 



GS98 AGS05 





l ini 


1 surf 


y 
* ini 


1 surf 


Solar param. 


OLi 


a 
Pi 




a 
Pi 


sll 


0.0606 


0.1290 


0.0603 


0.1393 


s33 


-0.0037 


-0.0047 


-0.0034 


-0.0044 


s34 


0.0079 


0.0098 


0.0075 


0.0095 


sl7 


0.0000 


0.0000 


0.0000 


0.0000 


sbe7e 


0.0000 


0.0000 


0.0000 


0.0000 


shep 


0.0000 


0.0000 


0.0000 


0.0000 


snl4p 


0.0010 


0.0011 


0.0006 


0.0007 


aue 


-0.1360 


-0.1954 


-0.1408 


-0.2002 


difu 


0.0271 


-0.0733 


0.0331 


-0.0733 


lumi 


0.3770 


0.3522 


0.4070 


0.3810 


C 


-0.0079 


-0.0073 


-0.0040 


-0.0020 


N 


-0.0017 


-0.0009 


-0.0004 


0.0010 





0.0083 


0.0183 


0.0140 


0.0285 


Ne 


0.0285 


0.0355 


0.0237 


0.0313 


Mg 


0.0277 


0.0317 


0.0314 


0.0366 


Si 


0.0593 


0.0582 


0.0682 


0.0696 


S 


0.0439 


0.0414 


0.0506 


0.0480 


Ar 


0.0123 


0.0112 


0.0075 


0.0070 


Fe 


0.0813 


0.0842 


0.0877 


0.0943 



Table 2. Initial solar helium abundance predictions from different standard solar models 



SSM 


^ini,0 


^surf.O 


v^ini 
I Q 


Villi 
I G 








K (l>s) = 0.836 


X ( i,s) = 0.903 


GS98 


0.2721 


0.2423 


0.2779 


0.2782 


AGS05 


0.2593 


0.2292 


0.2774 


0.2785 


AGSS09 


0.2617 


0.2314 


0.2777 


0.2787 


AGSS09ph 


0.2653 


0.2349 


0.2771 


0.2788 



